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Using recently developed methods for the evaluation of five- loop amplitudes in perturbative QCD, 
corrections of order oT s for the non-singlet part of the cross section for electron-positron annihilation 
into hadrons and for the decay rates of the Z-boson and the r-lepton into hadrons are evaluated. The 
new terms lead to a significant stabilization of the perturbative series, to a reduction of the theory 
uncertainly in the strong coupling constant a s , as extracted from these measurements, and to a small 
shift of the central value, moving the two central values closer together. The agreement between two 
values of a 3 measured at vastly different energies constitutes a striking test of asymptotic freedom. 
Combining the results from Z and r decays we find a s (Mz) = 0.1198 ± 0.0015 as one of the most 
precise and presently only result for the strong coupling constant in order a*. 

PACS numbers: 12.38.Bx, 13.35.Dx, 13.85.Lg 
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The strong coupling constant a s is one of the three 
fundamental gauge couplings constants of the Standard 
Model (SM) of particle physics. Its precise determination 
is one of the most important aims of particle physics. Ex- 
periments at different energies allow to test the predic- 
tions for its energy dependence based on the renormal- 
ization group equations, the comparison of the results 
obtained from different processes leads to critical tests 
of the theory and potentially to the discovery of physics 
beyond the Standard Model. Last but not least, the con- 
vergence of the three gauge coupling constants related 
by SU(3)xSU(2)xU(l) to a common value, after evolv- 
ing them to high energies, allows us to draw conclusions 
about the possibility of embedding the SM in the frame- 
work of a Grand Unified Theory. 

One of the most precise and theoretically safe determi- 
nation of a s is based on measurements of the cross section 
for electron-positron annihilation into hadrons. These 
have been performed in the low-energy region between 
2 GeV and 10 GeV and, in particular, at and around the 
Z resonance at 91.2 GeV. Conceptually closely related is 
the measurement of the semileptonic decay rate of the r- 
lepton, leading to a determination of a s at a scale below 
2 GeV. 

From the theoretical side, in the framework of pertur- 
bative QCD, these rates and cross sections are evaluated 
as inclusive rates into massless quarks and gluons [3, Q • 
(Power suppressed mass effects are well under control for 
e + e~-annihilation, both at low energies and around the 
Z resonance, and for r decays 0, |3,M @, S @|j an d the 
same applies to mixed QCD and electroweak corrections 

0E1). 



The ratio R(s) = a(e + e 



hadrons)/cr(e 



correlator of the electromagnetic current j^: 

R(s) = 127rlmn(-s-ie), (1) 
3Q 2 IT(Q 2 ) = i [ d 4 ^-(0|Tj AI (x)^(0)|0), (2) 



with Q 2 = —q 2 . It is also convenient to introduce the 
Adler function as 



D(Q 2 ) = -12tt 2 Q : 



^n(Q 2 ) - q (s+q2) 2 



R(s)ds 



R(s) = D(s) - ir 2 [3 2 {^al + (da + A-<k ft) o*} + 



We define the perturbative expansions 



D(Q 2 ) 
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i=0 



d^KQ 2 ), R(s) 
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nai(s), (3) 



where a s = a s /Tt and the normalization scale is set to 
/i 2 = Q 2 or to /j, 2 — s for the Euclidian and Minkowskian 
functions respectively. The results for generic values of 
ix can be easily recovered with standard RG techniques. 

Note that the first three terms of the perturbative se- 
ries for D and R coincide. Starting from ra, terms pro- 
portional 7T 2 arise which can be predicted from those 
of lower order. It has been speculated that these "7r 2 - 
terms" , also called "kinematical terms" , might constitute 
a major part of the full higher order corrections (see, e.g. 



ll|, [12[ and references therein); however, the validity of 



li + fi ) is expressed through the absorptive part of the 



this hypothesis can only be established by the full calcu- 
lation. Indeed, for the scalar correlator this assumption 
has been shown to fail 13 1. 

For the vector correlator the terms of order a 2 and 
have been evaluated nearly thirty and about fifteen years 
ago [3, EH, HH , respectively. The corrections are con- 
veniently classified according to their power of n/ , with 
ro/ denoting the number of light quarks. The a^n^ term 
is part of the "renormalon chain" , the evaluation of the 



2 



next term, of order a\r?p was a test case for the tech- 
niques used extensively in this paper and, furthermore, 
led to useful insights into the structure of the perturba- 
tive series already [17 1. 

The complete five-loop calculation requires the evalu- 
ation of about twenty thousand diagrams (we have used 
QGRAF [l8| for their automatic generation). Using "in- 
frared rearrangement" [l9| , the R* operation [2(1 and the 



prescriptions formulated in [2l| to algorithmically resolve 
the necessary combinatorics, it is possible to express the 
absorptive part of the five- loop diagrams in terms of four- 
loop massless propagator integrals. 

These integrals can be reduced to a sum of 28 mas- 
ter integrals with rational functions of the space-time di- 
mension D as coefficients. The latter ones were fully 
reconstructed after evaluating sufficiently many terms of 
the 1/D expansion [23j of their representation proposed 



in [22j. This direct and largely automatic procedure re- 
quired enormous compu ting resources and was performed 
using a parallel version [2J 



25|. 



of FORM 

In this paper we present the results for the so-called 
"non-singlet" diagrams. These are sufficient for a com- 
plete description of r-decays. For e + e _ annihilation 
through a virtual photon they correspond to the dom- 
inant terms proportional X^Qf- The singlet contribu- 
tions proportional (J2iQi) 2 arise for the first time in 
(D(a.g). They are known to be small, and will be eval- 
uated at a later point. Similar comments apply to the 
singlet contributions in Z decays. 

The analytic result for the five loop term in the 
Adler function is given by (we suppress the trivial fac- 
tor 3 f Q "f throughout) 



ck = n 3 f [- 

,2 r 1045381 
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5832 
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The knowledge of d 4 leads straightforwardly to R at 
order oq, for brevity given below in numerical form: 



R 



1 + a s + (1.9857 - 0.1152 n/) a 2 s 
(-6.63694-1.20013n/ 



0.005187^)^ 



(5) 



+ (-156.61 + 18. 77n/-0.7974n : 



0.0215 nl) a* 



It is also instructive to explicitly display the genuine five- 
loop contributions to d 4 (underlined in (|6I7[) ) and the 
"kinematical" terms originating from the analytic con- 
tinuation: 



r3 = 



18,2 - 24.9 + ( -4.22 + 
( -0.086 + 0.091) nj, 
135^8- 292.4+ ( :z 34 I 4 
( 1.88 - 2.67) n) + ( -0.010 



3.02) n s 
+ 53.2)n/ 



(0) 
(7) 



0.032) nl. 



Since it will presumably take a long time until the next 
term of the perturbative series will be evaluated, it is of 
interest to investigate the predictive power of various op- 
timization schemes empirically. Using the principles of 
"Fastest Apparent Convergence" (FAC) [HI or of "Min- 
imal Sensitivity" (PMS) [2jJ, which happen to coincide 



in this order, the central values of the predictions 11 



df ed (n f = 3,4,5) = 27+16,8 + 18,-8 + 44 

differ significantly from the exact result, 

dT ac \n f = 3,4,5) = 49.08, 27.39, 9.21. (8) 

However, within the error estimates [28j, predicted and 
exact values are in agreement. The picture changes, once 
these estimates are used to predict the coefficient r 4 . Al- 
though sizable cancellations between "dynamical" and 
"kinematical" terms are observed for the individual n/ 
coefficients in ([7]) the predictions for the final results are 
significantly closer (in relative sense) to the results of the 
exact calculation: 



prcd 



(n/ = 3,4,5) = 
t (n f = 3,4,5) = 



-129 + 16, -112 + 30,-97 + 44, 
-106.88, -92.898, -79.98. 



(This is in striking contrast to the case of the scalar cor- 
relator, where the predictions for the dynamical terms 
work well, but, as a consequence of the strong cancel- 
lations between dynamical and kinematical terms fail in 
the Minkowskian region [l3(.) 

Using FAC and the exact result for ch, the coefficients 
rf 5 and r 5 can be predicted (following [lj [5§| ) for rif — 
3, 4, 5, namely 



dl Iea (n f = 3,4,5) = 275, 152, 89, 



(9) 



rr d (n/ = 3,4,5) 



-505, -134, 168. (10) 



These terms may become of relevance for the In- 
ternational Linear Collider (LLC) running in the 
GIGA-Z mode with an anticipated precision of 
5a s = 0.0005 — 0.0007 0, and already today for the 
analysis of r-decays. 

From the combined analysis of data for er(e + e~ — > 
hadrons) in the region between 3 and 10 GeV a value 



a s {9 GeV) = 0.182 + 0.033 



(11) 



has been obtained recently [30( . The shift in a s from the 
inclusion of the a 4 s term amounts to <5a s (9 GeV) = 0.003 
and is thus irrelevant compared to the large experimental 
error. 

The situation is different for Z decays. The analysis 
of the electroweak working group [3l[ is based on eq. (J5[) 
with nt = 5, including term up to 0(oq) and leads to 



a s (M z ) 



NNLO 



0.1185 + 0.0026 



(12) 
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Since additional corrections (mixed QCD /electroweak or 
mass terms) are only weekly a s -dependent we may con- 
sider R(s = M§) as a pseudo-observable: 

R(s = M§) = 1.03904 ± 0.00087. (13) 

Including the a* term leads to a shift 5a s (Mz) — 0.0005 

a s (M z ) NNNLO = 0.1190 ± 0.0026 cxp , (14) 

The theory error may either be conservatively based on 
the shift produced by the last term (0.0005) or on the 
scale variation with fi/Mz = | -f- 3, leading to ±0.0002 
and can be neglected in both cases. 

Higher orders are of particular relevance in the low- 
energy region, for example in r decays. The correction 
from perturbative QCD to the ratio 



R T 



V+A 



F(t — >hadronS5 = o+7y T ) 



i\y U d\ 2 SEw{^ 



7$ 

S 



6' 



(15) 



EW 



JNP) 



is given by 



l+<5 



In the subsequent analysis we will use Sew = 1-0198 ± 
0.0006 and S' EW = 0.001 for the electroweak corrections, 
8 2 = (-4.4 ± 2.0) x 10~ 4 for light quark mass effects, 
Snp_= (— 4.8±1.7) x 10 -3 for the nonperturbative effects 
@,E3 and V ud = 0.97418 ± 0.00027 0. 

The perturbative quantity 5q can be evaluated in Fixed 
Order perturbation theory or with "Contour Improve- 
ment" as proposed in 3J, l35j 



rFO _ 

°o — 

?C7 



5.202 at + 26.366 < 



Sfc 1 = 1.364 a s 



2.54 at + 9.71 at 



127.079 a 4 , (17) 
-(- 64.29 a 4 . (18) 



(To obtain the a s -dependent coefficients in eq. (TT8| we 
follow [33, IM 0, H3| and use a s (M T ) = 0.334 as refer- 
ence value.) For the subsequent analysis we will use as 
starting point <5g Xp = 0.1998±0.0043 exp as obtained from 
[36| and R T y+A = 3.471 ± 0.011 which in turn is based 
on the " universality- improved" electronic branching ratio 
B e = (17.818±0.032)% and the world average of the ratio 
of strange hadronic and electronic widths 0.1686±0.0047. 
The new values of a s {M T ) in dependence on the choice of 
di (with the previous estimate and the new exact result) 
are summarized in Table I. 

As stated above the theory error for a s from Z decays 
is small compared to the experimental uncertainties. The 
situation is more problematic for t decays and to some 
extent the theory error remains to be matter of choice. 



As anticipated in [28j it decreases significantly, once a* 
terms are included. However, the difference between the 
two methods stabilizes (this was checked in [28] by adding 



TABLE I: Results for a s (M T ) for different values of d A . The 
first line displays the a 3 s results (with the a 4 terms set to zero 
in eqs. (|17I18|| ). The second line uses the the previously pre- 
dicted value for d^, the last one uses the exact result Q. The 
first error is the experimental one; the second (theoretical) 
uncertainty in the value of a s corresponds to changing the 
renormalization scale /x as follows m7 M t = 0.4 - 2. 







af°(M T ) 


qP(m t ) 






0.337 ± 0.004 ± 0.03 


0.354 ± 0.006 ± 0.02 




= 25 


0.325 ± 0.004 ± 0.02 


0.347 ± 0.006 ± 0.009 


(1.4, 


= 49.08 


0.322 ±0.004 ±0.02 


0.342 ± 0.005 ± 0.01 



an estimate for the term) and must be considered as 
irreducible uncertainty. Given the input specified above 
we obtain as our final result 



a s {M T ) = 0.332 ± 0.005 cxp ± 0.015 



5th- 



(19) 



For the central value we take the mean value of FO and 
CI. For the theory error we take half of the difference be- 
tween two methods (0.01) plus (module of) the estimated 
correction from a]l term (-0.005), the latter being based 
on d 5 = 275 (see eq. ©). 

Applying four-loop running and matching [38l . 39L 40l . 
|4~D] to HH) we arrive at 

a s (M z ) = 0.1202 ±0.0006 oxp ± 0.0018th ± 0.0003 cvo i , 
= 0.1202 ±0.0019. (20) 

Here the evolution error receives contributions from the 
uncertainties in the c-quark mass (0.00003, m c (m c ) = 
1.286(13) GeV 0) and the 6-quark mass (0.00001, 
m b (m b ) = 4.164(25) GeV [42]), the matching scale 
(0.0001, /x varied between 0.7 m q (m q ) and 3.0 m 9 (m g )), 
the four-loop truncation in the matching expansion 
(0.0001) and the four-loop truncation in the RGE equa- 
tion (0.0003). (For the last two errors the size of the shift 
due the highest known perturbative term was treated as 
systematic uncertainty.) The errors are added in quadra- 
ture. 

Summary: The exact result for the a 4 term in the 
Adler function allows to extract the strong coupling con- 
stant from Z and r decays with high precision. Including 
the exact a 4 leads to small shifts of the central values 
and to a significant reduction of the theory uncertainty. 
Note that the shifts in a s (Mz) from Z- and r-decays, 
are opposite in sign and move the values in the proper 
direction, decreasing, thus, the current slight mismatch 
between two independent determinations of a s . 

The final results 

a s (M z )\z = 0.1190 ±0.0026, 
a s (M z )\ r = 0.1202 ±0.0019 



4 



from these two observables, although based on measure- 
ments of vastly different energy scales, are in remarkable 
agreement. This constitutes a striking test of asymptotic 
freedom in QCD. The two values can be combined to 

a s (M z ) = 0.1198 ±0.0015. 

This is one of the most precise and presently only result 
in order a*. 
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